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ABSTRACT 

We present a comparison of the SCUBA Half Degree Extragalactic Survey (SHADES) at 450 //m, 850yum and llOOyum with deep 
guaranteed time 15 /im AKARI FU-HYU survey data and Spitzer guaranteed time data at 3.6 - 24 yum in the Lockman Hole East. The 
AKARI data was analysed using bespoke software based in part on the drizzling and minimum- variance matched filtering developed 
for SHADES, and was cross-calibrated against Infrared Space Observatory (ISO) fluxes. Our stacking analyses find AKARI 15 fim 
galaxies with ; 200 /iJy contribute > 10% of the 450 yum background, but only < 4% of the llOOyum background, suggesting that 
different populations contribute at mm- wavelengths. We confirm our earlier result that the ultra-deep 450 yum SCUBA-2 Cosmology 
' Survey will be dominated by populations already detected by AKARI and Spitzer mid-infrared surveys. The superb mid-infrared 

' wavelength coverage afforded by combining Spitzer and AKARI photometry is an excellent diagnostic of AGN contributions, and we 

' find that (23 - 52)% of submm-selected galaxies have AGN bolometric fractions /agn > 0.3. 
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j_j 1 . Introduction adjusting the physical parameters used to characterised feedback 

processes (e.g. Granato et al. 2006) and/or with adjustments to 



the initial mass function (e.g. Baugh et al. 2005, Lacey et al. 
2008). 



The pioneering submm-wave observations of the Hubble Deep 
Field North (Hughes et al. 1998, Barger et al. 1998) and galaxy 
cluster gravitational lenses (Small et al. 1997) revealed the Indeed the existence of strong relationships between super- 
presence of completely unanticipated populations of submm- massive black hole masses in nearby galaxies and the proper- 
luminous galaxies with fluxes of several mJy at 850//m. This ties of the host bulges or spheroids (e.g. the Magorrian rela- 
discovery was one of the first in a series that led to the devel- tion, Magorrian et al. 1998; see also e.g. Ferrarese & Merritt 
opment and ultimately wide acceptance of "downsizing" phe- 2000) all but require a strong link between high-redshift black 
nomenological models of galaxy evolution, that massive galax- hole accretion and the generation of the stellar mass of the host, 
ies formed the bulk of their stars at high redshifts in large star- This is reinforced by the observations that many or perhaps all 
burst events, with the opposite behaviour found in lower-mass quasars show evidence of star formation (e.g. Lehnertet al. 1992, 
systems. This is contrary to naive expectations from hierarchical Hughes et al. 1997, Aretxaga et al. 1998, Brotherton et al. 1999, 
models, but nevertheless it was found these observations could Lutz et al. 2008, Mullaney et al. 2009, Shi et al. 2009, Veilleux et 
be made consistent with semi-analytic hierarchical models by al. 2009, Serjeant & Hatziminaoglou 2009). While submm and 
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mm-wave direct detections are only made in a minority of high- 
redshift quasars (e.g. Carilli et al. 2001, Omont et al. 2001, Isaak 
et al. 2002, Bertoldi & Cox 2002, Omont et al. 2003, Priddey et 
al. 2003a,b, Willlot et al. 2003, Bertoldi et al. 2003), stacking the 
far-infrared to mm-wave non-detections yielded significant de- 
tections and demonstrated clearly that all quasars are on average 
ultraluminous starbursts (Serjeant & Hatziminaoglou 2009). If 
submm-selected galaxies are the violent starbursts that mark the 
sites of the progenitors of present-day giant ellipticals, then the 
presence and nature of AGN in submm-selected galaxies may 
track the feedback processes regulating the stellar mass assem- 
bly. 

Are the bolometric luminosities of submm galaxies domi- 
nated by the release of gravitational binding energy around a 
central supermassive black hole, or are they dominated by the re- 
lease of nuclear binding energy in starbursts? Similar questions 
were asked of local ultraluminous infrared galaxies in the 1990s 
(see e.g. Genzel & Cezarsky 2000); now the issues of feed- 
back and the Magorrian relation add an additional perspective 
and motivation to the higher-redshift equivalent questions for 
submm-selected galaxies. Deep 2 - lOkeV imaging of submm- 
selected galaxies has found that around three-quarters show evi- 
dence of an active nucleus (Alexander et al. 2005), the majority 
of which are heavily obscured. Nevertheless, the X-ray to far- 
infrared luminosity ratios suggest that star formation rather than 
active nuclei typically dominate the bolometric output, assuming 
the underlying spectral energy distributions of the active nuclei 
resemble those of unobscured quasars. In order to eliminate the 
possibility that Compton-thick active nuclei in submm-selected 
galaxies have different underlying spectral energy distributions 
than the comparator populations, one should sample the spectra 
around the expected peak of the dust torus output, in the mid- 
infrared. 

The advent of the Spitzer space telescope led to the dis- 
covery of new dust-shrouded galaxy populations with starbursts 
and/or active nuclei (e.g. Martmez-Sansigre et al. 2005, Dey 
et al. 2008). The mid-IR where the dust torus luminous out- 
put peaks is the ideal place to test for AGN; we expect star- 
forming galaxies to have strong PAH and silicate features, while 
AGN dust tori have much more featureless spectra, with the 
possible sole exception of a lOjum rest-frame silicate absorp- 
tion trough. We found in Negrello et al. 2009 that with suf- 
ficiently comprehensive broad-band photometry in the mid-IR, 
starburst galaxies can have photometric redshift determinations 
from the redshifting of PAH features. Furthermore, if the bolo- 
metric AGN fraction is more than 30%, the mid-IR SEDs are 
sufficiently close to featureless power-laws in broad-band pho- 
tometry that photometric redshifts are impossible with mid-IR 
data alone. However in these cases, the power-law SEDs are 
a useful signature of bolometrically-significant AGN regard- 
less of redshift. In this case, the longer wavelength range ob- 
tained by combining AKARJ and Spitzer offers significant ad- 
vantages over the insights available (albeit useful ones) from 
broad-band Spitzer photometry alone (e.g. Yun et al. 2008). A 
recent mid-infrared spectroscopic study of 13 submm-selected 
galaxies (Pope et al. 2008) found only 2/13 with such power-law 
mid-infrared spectral energy distributions; similarly, Menendez- 
Delmestre et al. (2009) found 4/23 submm-selected galaxies 
with continuum-dominated mid-infrared spectra. Clearly there 
is a need for larger samples of submm-selected galaxies with 
constrained AGN bolometric fractions. 

Submm galaxies are also significant contributors to the 
submm-wave extragalactic background light (EBL). The EBL 
can be expressed as an integral over the comoving volume- 



averaged luminosity density (with an additional (1 +z) factor, e.g. 
Peacock 1999), so at any redshift z the galaxies that dominate the 
far-IR to mm-wave EBL are also the ones that dominate the star 
formation rate at that z- But resolved submm galaxies only con- 
tribute a few tens of percent at most (Hughes et al. 1998) to the 
EBL at those wavelengths. To constrain the populations that con- 
tribute the rest one needs to use stacking analyses. Several pop- 
ulations of galaxies have been shown to contribute significant 
minorities of the submm EBL, such as EROs and HEROs, BzK 
galaxies, (e.g. Coppin et al. 2004, Takagi et al. 2007, Greve et 
al. 2009), and Lyman-break galaxies (e.g. Peacock et al. 2000). 
In Serjeant et al. 2008 we found that the 24 jum-selected galaxies 
contribute the bulk of the 450 /vm background, but that at 850 //m 
there are unknown new populations that contribute the bulk of 
that background. 

In this paper we will constrain the AGN bolometric frac- 
tions of a sample of submm-selected galaxies in the Lockman 
Hole East, taken from the SCUBA Half Degree Extragalactic 
Survey (SHADES), using mid-infrared photometry from the 
AKARI and Spitzer space telescopes. We also constrain the con- 
tributions the mid-infrared populations make to the submm and 
mm-wave extragalactic background light using stacking analy- 
ses. SHADES papers I and II (Mortier et al. 2005 and Coppin 
et al. 2006) present the SHADES survey design, data analy- 
sis and source counts. Paper III (Ivison et al. 2007) gave radio 
and mid-infrared counterparts of SHADES galaxies. SHADES 
paper IV (Aretxaga et al. 2007) presents photometric redshift 
constraints of submm galaxies derived from far-infrared and ra- 
dio data alone. SHADES paper V (Takagi et al. 2007) con- 
strains the submm properties of near-infrared-selected galax- 
ies. Paper VI (Coppin et al. 2007) gave 350;um photometry 
for SHADES galaxies, while papers VII and VIII (Dye et al. 
2008 and Clements et al. 2008 respectively) analysed the mul- 
tiwavelength spectral energy distributions based on data avail- 
able at the time. SHADES paper IX (Serjeant et al. 2008) made 
submm stacking analyses of ISO and Spitzer-selected galaxies, 
and in this paper we will extend this analysis to include the 
AKARI catalogues described in the section|2] Austermann et al. 
(2009) extended the SHADES survey to mm-wavelengths and 
substantially increased the areal coverage. The comparison of 
the 850 /im and 11 00 yum data, including a critical comparison 
of the 1 100 /im data from different instruments and a systematic 
search for populations of submm-wave drop-outs, will be given 
in Negrello et al. (2009, in preparation). In this paper we will 
use the AzTEC 1 100 im\ data, rather than the 1 100 fim data from 
Coppin et al. 2008 which has not been corrected for Eddington 
bias (Eddington 1913). We will also use the multi-wavelength 
data from SHADES papers III, VI, VII and VIII (Ivison et al. 
2007, Coppin et al. 2007, Dye et al. 2008 and Clements et al. 
2008 respectively) who presented spectral energy distribution 
(SED) fits of SHADES galaxies. 

2. Data acquisition 

2.1. Submm and mm-wave observations 

The SCUBA Half Degree Extragalactic Survey (SHADES) was 
the largest project on the James Clerk Maxwell Telescope from 
2003-2005. During the instrument lifetime of SCUBA the sur- 
vey mapped a total of approximately a quarter of a square de- 
gree at 850 /im to a typical 3.5 cr depth of 8 mJy, in the Lockman 
Hole East and the Subaru-XMM Deep Field (Mortier et al. 
2005). The originally-planned areal coverage and more was 
completed by the SHADES team with the AzTEC instrument 
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at llOOyum (Austermann et al. 2009) to typical 3.5<t depths in 
the Lockman Hole East of 3.1 - 4.6 mJy. Both SHADES fields 
benefit from enormous multi-wavelength wide-field campaigns, 
and the Lockman Hole East in particular has some of the deepest 
Spitzer Space Telescope mid-infrared mapping of any contigu- 
ous field over hundreds of square arcminutes. More details on 
the SHADES project, including the survey goals, methodology 
and submm data analysis, can be found in SHADES papers I 
and II (Mortier et al. 2005, Coppin et al. 2006). Details on the 
SHADES mm-wave data analysis are given in Austermann et al. 
(2009). We use the Austermann et al. (2009) mm-wave measure- 
ments in preference to Coppin et al. (2008), since the former but 
not the latter have been corrected for Eddington bias. 

2.2. Mid-infrared observations 

The Lockman Hole East was targetted by the Spitzer space tele- 
scope in guaranteed time. The field was mapped with the IRAC 
and MIPS instruments (Fazio et al. 2004, Rieke et al. 2004) and 
sources extracted to depths of 4.47//Jy at 3.6 ;um (3cr), 4.54//Jy 
at 4.5 yum (3cr), 20.9yuJy at 5.8 fim (3cr), 12.5yuJy at 8 yum (3.2cr) 
and 38/iJy at 24 jum (4cr). For more details of this catalogue, see 
Serjeant et al. 2008 and references therein. The Lockman Hole 
East was also mapped at 15;um with the CAM instrument on 
the Infrared Space Observatory (ISO). More details of the ISO 
survey can be found in Elbaz et al. 1999 and Rodighiero et al. 
2004. 

The Lockman Hole East was also the subject of a 10' x 30' 
survey in the L15 15 /im filter by Infrared Camera (IRC) of 
AKARI, as part of the FU-HYU mission programme (Pearson et 
al. 2009), to a typical Icr noise level of 2Qfi]y. This programme 
targetted well-studied Spitzer fields to capitalize on the com- 
prehensive multi-wavelength data sets available in these fields. 
Some FU-HYU observations in the GOODS-N field have been 
presented in Negrello et al. (2009). 

A pipeline analysis of the IRC is available in the Image 
Reduction and Analysis Facility (IRAF). Despite being an ex- 
cellent general-purpose analysis, we found that the processing 
was not optimal for our purposes of minimum-variance point 
source extraction. The standard pipeline locates and flags bad 
pixels, then interpolates the image to correct for field distortions. 
However this interpolation stage creates Moire effects, correlates 
the noise between pixels and in some cases spreads bad individ- 
ual pixel readouts over several pixels in the interpolated image. 
Furthermore, some of the AKARI observations had been mistak- 
enly taken without the requested dithering intended to remove 
bad pixels, so our only option was to use the spacecraft jitter, 
which necessarily requires a bespoke analysis. 

We had previously developed mapping and minimum- 
variance point source extraction tools in the Interactive Data 
Language (IDL) for the SHADES survey, so we opted to use 
a similar approach here. We stopped the IRC pipeline after the 
flat fielding, imported the data into IDL, then determined an ap- 
proximate pointing solution using catalogues of ISO and Spitzer 
sources in the field. The relative jitter was found by measur- 
ing centroids of bright sources in each data frame, and hence 
we could determine a World Coordinate System solution for 
each frame. Bad pixels were identified and masked using the 
same pixel histogram fitting procedure employed in the SCUBA 
Hubble Deep Field North and SHADES observations (Serjeant 
et al. 2003a, Mortier et al. 2005). Several frames were affected 
by Earthshine, so we performed a polynomial fit to each row, 
then each column, then identified and corrected a small number 
of pixels which reported consistently low counts regardless of 



target. This latter effect may indicate the need to improve the flat 
field determinations, but the offsets induced by these corrections 
are much smaller than our photometric errors. Following the 
best practice in SCUBA observations we determined the noise in 
each detector pixel using Gaussian fits to readout histograms. We 
adopted the field distortion corrections as coded in the IRAF IRC 
pipeline, but from a comparison of our initial AKARI catalogues 
against Spitzer observations we found evidence for further field 
distortions. We fit to these using a second-order polynomial and 
incorporated the correction into our total field distortion solu- 
tion. 

To mosaic the individual frames onto a master image we used 
the "zerofootprint" drizzling technique developed originally for 
SCUBA (Serjeant et al. 2003a): detector pixels are mapped onto 
much finer sky pixels, accounting for field distortion and jitter, 
with the flux from a detector pixel deposited in a single sky pixel. 
Where a sky pixel has multiple observations, fluxes are com- 
bined using noise-weighted coadds. The final coadded map rep- 
resents a detector pixel's view of the sky, in the sense that at ev- 
ery point it reports the noise-weighted mean flux of all detector 
pixels centred exactly on that point. Notwithstanding cross-talk 
between detector pixels (which we neglect), the final coadded 
map also has statistically independent pixels, i.e. without covari- 
ances. 

We performed a further deglitching stage on the coadded 
map. The minimum- variance estimator for the point source flux 
F at any point in the map can be expressed a noise-weighted 
point spread function convolution, as used in submm surveys 
(Serjeant et al. 2003a, Mortier et al. 2005): 



F = 



(IW) 8) P 



(1) 



where / is the coadded image, W is the weight image (the recip- 
rocal of the noise squared), and P is the point spread function. 
Propagating errors on equation [T] yields an expression for AF, 
the error on F: 



(AFf = 



1 



(2) 



For the point spread function we assumed a Gaussian with a full 
width half maximum of 5.47". 

Our point source detection map is the F/AF image, in which 
we found objects using a connected pixels algorithm. Objects 
were identified as discrete regions with map values greater than 
a given threshold. We used a series of thresholds of 3 - lOcr in 
steps of Icr, because objects which are blended at a low thresh- 
old may be distinct objects at a brighter threshold (e.g. Mortier 
et al. 2005), yielding a catalogue of 622 distinct objects. The 
IRC pipeline flux calibration is given only for photometry taken 
with the IRC pipeline outputs, which is necessarily systemati- 
cally different to our pipeline outputs, so we cross-matched our 
catalogue with the previous ISO 15 /urn sources from Rodighiero 
et al. (2004) and based our calibration on the cross-matchings so 
the average flux ratio is unity, shown in Figs.[T] The final AKARI 
15 jum sources are shown in Figs.|2]and[3] 

3. Methods 

For the stacking analyses of submm galaxies we adopted the 
same methodologies of Serjeant et al. 2008. In brief, we first 
subtract the catalogued objects from the 850jum and 1100//m 
maps (there are no secure 450yum point source detections). We 
used Gaussian point spread functions for the submm data and the 
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Fig. 1. Comparison of 15 fim fluxes obtained in the Lockman 
Hole East with the Infrared Space Observatory (Elbaz et al. 
1999, Rodighiero et al. 2004) with those obtained with our 
AKARI survey (Pearson et al. 2009). Note that the mean flux cal- 
ibration of AKARI has been determined by the cross-matching 
with the ISO survey (implicitly using the ISO flux calibration) 
and is not determined independently in this study. 
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Fig. 3. The SHADES AzTEC llOOyum image of the Lockman 
Hole East. The AKARI 15 fim sources are overplotted in yellow 
(rectangular coverage) and the ISO 15 fim sources are overplot- 
ted in blue (approximately square coverage). 
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Fig. 2. The SHADES SCUBA 850 jum signal-to-noise image of 
the Lockman Hole East. The AKARI 15/im sources are over- 
plotted in yellow (rectangular coverage) and the ISO 15/vm 
sources are overplotted in blue (approximately square coverage). 



point spread function in Austermann et al. (2009) for the AzTEC 
data. We refer to these maps as the residual maps. The compar- 
ison of the histograms of mm-wave fluxes in the map and at the 
positions of 15//m-selected galaxies yielded a possibly strong 
signal (~ 8cr) but which appeared to be dominated by a non- 



Gaussian shape to the latter histogram (Fig. |4]i. This could be 
suggestive of weak but coherent features in parts of the map, but 
which would necessarily have to have absolute values less than 
Icr per beam, making investigation difficult. Alternatively, this 
could be a consequence of having a relatively small number of 
AzTEC beams in the 15;um survey coverage. Without passing 
judgement on these issues, we have opted to take a very conser- 
vative approach and adopted an additional beamswitching sky 
subtraction, by subtracting from the AzTEC residual map the 
same map offset by 100 pixels East, or 300". We also derived a 
noise map for this beamswitched residual map by measuring the 
standard deviation in a 45" radius circular region at every point 
in the map. 

There is some disagreement in the literature over the best 
practice in stacking methodology Some (such as Peacock et al. 
2000 and Serjeant et al. 2004, 2008) have opted to subtract 
known point sources, while others (such as Wang et al. 2007, 
Scott et al. 2008, Marsden et al. 2009) opt to use the unsubtracted 
map. In escense, it depends on what one wants to measure. If 
one requires the total intensity from all galaxies, one should not 
subtract any flux (or the subtracted population should be added 
afterwards). If however one requires the mean flux from a par- 
ticular population, as distinct from all galaxies as a whole, then 
one needs to clean the map to avoid any biases in the mean. 

In Serjeant et al. (2008) we found that the submm-selected 
galaxies had very different mid-IR/submm flux ratios to the 
galaxies selected at 24yum. This could be due to bimodality, i.e. 
two populations, or we could be sampling two ends of a contin- 
uum; it is difficult to distinguish these possibilities from stacks 
alone. The possibility that a few 15 /im-selected galaxies are (rel- 
ative to the rest) pathological would suggest the use of point- 
source-subtracted maps, which is the approach we have adopted. 

We take unweighted mean averages of the submm fluxes in 
the residual maps (or for the AzTEC observations, mm-wave 
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fluxes) at the positions of mid-infrared-selected galaxies. By the 
Central Limit Theorem the unweighted mean flux is Gaussian 
distributed with a dispersion of cr/ Va^ where cr is the standard 
deviation of the sample of fluxes and is the number of fluxes 
being averaged. We also use Kolmogorov-Smirnoff' tests to com- 
pare the histogram of fluxes at the positions of mid-infrared- 
selected galaxies with the histogram of the submm/mm-wave 
fluxes in the submm/mm-wave map as a whole, excluding re- 
gions with no mid-infrared data. This latter test inherently in- 
corporates a control comparison. We include the 15//m survey 
from ISO (Elbaz et al. 1999, Rodighiero et al. 2004); although 
this survey is not as deep as our AKARI mapping, it increases the 
areal coverage of the stacking analysis (Figs.|2]and[3]l and there- 
fore also the number of submm and mm-wave beams contribut- 
ing, which improve the stacked signal (Serjeant et al. 2008). 
Unweighted means give cosmetically poor results for stacked 
postage stamps due to edge effects, so we used noise-weighted 
means for the purposes of making stacked postage stamps only. 

For the characterisation of the spectral energy distributions 
of submm-selected galaxies, we assume a minimum photomet- 
ric error of ten percent. We use an ensemble of starburst model 
outputs from Takagi et al. (2003, 2004) and an ensemble of dust 
torus model outputs from Efstathiou & Rowan-Robinson (1995), 
following the methodology of Negrello et al. 2009. We make lin- 
ear superimpositions of these model spectral energy distributions 
neglecting the small effects of energy transfer between the dust 
torus and the starburst dust components. The key physical pa- 
rameters in these models are: the torus opening angle 6 (which 
we set to 45°), the ratio of outer and inner torus radii (which we 
set to 20), and viewing angle flyiew; starburst compactness pa- 
rameter 0, age T and dust composition (MiUcy Way, SMC or 
LMC). A star formation timescale of 100 Myr was adopted for 
the Takagi et al. models. A wider parameter space of AGN mod- 
els is not justified by the contraints available from our data. We 
leave the redshift as a further free parameter even where spec- 
troscopic redshifts have been published, to avoid dependence on 
the often-uncertain discussion of identifications. The total num- 
ber of fitted parameters (including normalisation, but not count- 
ing dust composition) is three, and two more with the addition 
of AGN, plus a further fitting parameter for redshift. 

4. Results 

4. 1 . Stacking analysis results 

Fig. |5] shows noise-weighted stacked 450/im, 850 ;um and 
1 100 ijm postage stamps of the AKARI and ISO 15 yum-selected 
galaxies. This is a significant improvement over the detections in 
Serjeant et al. (2008), and for the first time extends the stacked 
signal of mid-infrared-selected galaxies to mm- wavelengths. 
The unweighted mean submm-wave and mm-wave fluxes of 
the 15jum-selected galaxies are (54go^m) - (3.7 + 1.1) mJy, 
<5850^m> = (0.206 ± 0.084) mJy and (5 iioo^m) = (0.148 ± 
0.052) mJy. The 450 : llOOyum flux ratio is consistent with a 
grey -body index of [i — 1 .5 on the Rayleigh-Jeans tail, though 
the 850 yum stacked flux is marginally discrepant. In the source 
count model of Pearson (2005) the predicted median redshift of 
our sample is z ~ 1, with of the order 10% at z > 2. 

Fig.|4]shows the histograms of submm and mm-wave fluxes 
at the positions of the 15 /vm-selected galaxies, compared to the 
histograms of the maps as a whole were 15 jjm data is available. 
The probabilities that the 15 yum-selected galaxies are drawn ran- 
domly from the submm maps are ~ lO"-' at 450 jum, 0.046 at 
850;um and 0.016 at 1100//m, i.e. a significant stacking detec- 



tion. These probabilities in all three cases are due to slight asym- 
metries or positive shifts in the submm/mm-wave flux distribu- 
tions. 

4.2. SED fitting results 

We found centroids in the 15 jum AKARI point source detection 
map at the positions of SHADES Lockman galaxies reported in 
Dye et al. (2008). Where no centroid solution could be obtained, 
we took the value of the point source flux F and its associated er- 
ror AF at the SHADES counterpart position (equations[T]and|2]l. 
Tables[T]and|2]give a compilation of available photometry of the 
SHADES galaxies from Dye et al. (2008), Ivison et al. (2007) 
and Coppin et al. (2007), and incorporating our new 15 fim pho- 
tometry. 

Following Negrello et al. (2009), we found the minimum 
solution for the available photometry assuming a mix of starburst 
and active nucleus, allowing the relative bolometric fractions of 
the components (/sb and /agn = 1 - /sb respectively) to vary. 
Table |3] gives the best-fit parameters of these two-component 
fits, excluding those for which no fit could be made, but in- 
cluding for completeness those for which the fitting is under- 
constrained. Figs. l6]to[T0l show the best-fit SEDs, as well as the 
best-fit starburst-only model for each SHADES galaxy. In many 
cases, adding the active nucleus is a requirement for reproducing 
the mid-infrared data. In several other cases however, no accept- 
able could be found, because the submm and mm-wave data 
exceeded the range given by the models. In these cases, a fur- 
ther cold cirrus component may contribute significantly to the 
longest wavelengths. 

Several galaxies are worth discussing in detail. The galaxies 
with an excess at submm/mm-wavelengths are LOCK850.009, 
Oil, 013, 031, 037B, 040, 043/043B, 053, 070, and 071. In 
addition, Lock850.076 has peculiar submm/mm-wave colours, 
perhaps due to incorrect flux deboosting in either the submm 
or mm-wave data. The 15;um flux in Lock850.003B may 
be contaminated by a nearby object. The SED fits in some 
cases suggest particular identifications, e.g. LOCK850.010B 
is preferred over LOCK850.010 (not plotted since no accept- 
able fit could be made to the data), as is LOCK850.015 
over LOCK850.015B, LOCK850.035 over LOCK850.035B, 
and LOCK850.073 over LOCK850.073B. The identifica- 
tions for LOCK850.004, LOCK850.009, LOCK850.043 and 
LOCK850.77 remain ambiguous. 

5. Discussion and conclusions 

Following the methodology of Serjeant et al. (2008), we can es- 
timate the contribution that clustering of the mid-infrared popu- 
lation makes to the submm and mm-wave stacks, by integrating 
the correlation function w{6) with respect to angle 0: 

JfOO 
w(e)Bi6)2ned6 (3) 


where B is the submm or mm-wave beam and / is the total back- 
ground intensity. We argued in Serjeant et al. (2008) that the 
clustered contributions at 450 fim and 850 fim can be neglected. 
At 1100 //m equation 13 gives S = O.12(0o/arcsec)*' '* for an as- 
sumed w(ff) = (OlOoY^'^. The value of 9q is not well-determined 
for the mid-infrared population at these flux densities (typically 
; 200A(Jy), but Oliver et al. (2004) measured 9q = 1.24" for the 
faintest 3.6 //m-selected populations (note that nearly all 15 jim- 
selected galaxies are detectable at 3.6 //m with Spitzer in even 
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Table 1. Optical and near-infrared AB magnitude photometry of SHADES galaxies in the Lockman Hole East. Blank spaced 
indicated non-detections or lack of observations. The data are taken from Dye et al. 2008. 
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Table 2. Mid-infrared to radio photometry of SHADES galaxies in the Lockman Hole East. The 15 jum data are presented for the 
first time here, and the remaining data are taken from Coppin et al. 2006, Austermann et al. 2009, Ivison et al. 2007, Coppin et al. 
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Table 3. Results of the SED fitting, with the parameters defined in the text. All parameters are dimensionless except the age which 
is in Myr. In the final column, "C" notes possibly contaminated flux at 15 fim, "D" notes a depreciated identification, "I" insufficient 
data to determine the AGN bolometric fraction, "U" an uncertain identification, and "X" a far-infrared excess. 
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Fig. 4. Histograms of the submm and mm-wave fluxes in the regions with 15;um data (unhatched data, right-hand ordinates) 
compared to the submm and mm-wave flux measurements at the positions of 15 /im point sources (hatched data, left-hand ordinates). 
The significance values of Kolmogorov-Smirnov tests are quoted in the figures. As noted in the text, the 1 100 /jm map has had an 
additional sky subtraction from beamswitching. 



450jjm stack of ISpm galaxies 850pm stack of 15Mm galaxies 110DiJm stack of 15ijm galaxies 



o 
o 




Fig. 5. Stacked mm-wave and submm signal-to-noise postage stamps of 15 fim sources from AKARI and ISO, at (from left to right) 
450/zm, 850/im and 1100 //m. The images are scaled from -2.81cr to 3.73(T (450jum), -2.73cr to 2.83cr (850//m) and -3.01cr to 
2.99cr(1100yum). 
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the deepest 15jum exposures, e.g. Hopwood et al. 2009). To 
avoid introducing model-dependences in our quoted results, we 
opt not to correct for clustering in the stacked flux numbers 
quoted here. Nevertheless, the total background contributed by 
the 15 /im population in this paper is considerably less than the 
integrated background, so the clustered contribution estimated 
by this methodology is a very conservative upper limit. 

The observed 1100 yum stacked signal of (5iio()//m) = 
(0.148 + 0.052) mJy at a 15yum source combined with our ob- 
served surface density of 4950 ± 170 15 /im galaxies per square 
degree (Poisson errors) corresponds to a 1100//m background 
contribution of (0.73 ± 0.26) Jy deg"^. The estimated total back- 
ground at 1100/im is around ISJydeg"^ (e.g. Gispert et al. 
2000). Given the contribution from clustering, we can set an up- 
per limit of 4% to the contributions from these 15//m-selected 
galaxies to the 1100/im extragalactic background light. Our 
15 //m-selected sample is not as deep as the 24/im-selected pop- 
ulation used by Serjeant et al. (2008), but following the same 
methodology as we have applied at 1 100 /im, we estimate back- 
ground contributions of (1.02 ± 0.42)Jydeg"^ at 850 /im and 
18.3 ±5.4 Jy deg"^ and 450/^m, corresponding to (3 + 1)% of the 
850 //m background and (11 + 3)% of the 450 /(m background. 
This is consistent with our results in Serjeant et al. (2008), im- 
plying the populations which dominate the ultra-deep confusion- 
limited maps expected with the SCUBA-2 instrument will over- 
lap substantially with populations already detected with AKARI 
and Spitzer. 

The small background contribution at 1100//m is in stark 
contrast to our earlier results for the submm backgrounds 
(Serjeant et al. 2008) in which we found the majority of the 
450 yum background is attributable to the 24 /im-selected popula- 
tion, and to claims from the Balloon-Borne Large Area Submm 
Telescope (BLAST) data that the mid-infrared population is re- 
sponsible for all the 250 - 500 /im background (Devlin et al. 
2009, Marsden et al. 2009). Taken in combination with the re- 
sults in Serjeant et al. (2008), our results suggest strongly that 
other populations undetected by AKARI or Spitzer must domi- 
nate the 850 jim and mm-wave background. 

Not all our sample has sufficient mid-infrared data to dis- 
tinguish starburst and AGN mid-infrared contributions, but ta- 
ble [3] gives the results where available. In total 41 SHADES 
galaxies have sufficient data for this constraint, not counting 
multiple identifications. Of these, we find that our models can- 
not reproduce the far-infrared luminosity in 10 unambiguous 
identifications (LOCK850.013, 031, 037B, 038, 040, 053, 060, 
070, 071, 075) and in two further ambiguous identifications 
(LOCK850.009, 043). One possibility is that the spectral energy 
distributions have an additional cool cirrus component not ac- 
counted for in the models, which would further reduce an AGN 
bolometric contribution. Another possibility is ffiat ffie identi- 
fications are wrong, and the submm emission comes from ob- 
jects with larger far-infrared to optical/near-infrared luminosity 
ratios. Of the remaining (apparently) unambiguously-identified 
objects, 12 have AGN fractions > 0.3 (LOCK850.001, 003, 016, 
022, 027, 028, 041, 047, 052, 067, 079, 087), while 14 have star- 
burst fractions > 0.7 (LOCK850.010B, 014, 015, 024, 030, 033, 
034, 048, 066, 073, 076, 078, 083, 100). For LOCK850.004 and 
LOCK850.077 the two candidate identifications result in one 
starburst (/agn < 0.3) and one AGN (/agn > 0.3) interpreta- 
tion for each galaxy. Finally two further cases merit individual 
attention. In LOCK850.060, the addition of the AGN compo- 
nent created a worse fit in the submm; we opt to attribute greater 
weight to this long-wavelength data point and include this galaxy 
among the starburst. In LOCK850.076 there is additional mm- 



wave flux, which we attribute to errors in the deboosting ffie 
submm or mm-wave data in this individual object; this galaxy 
is also counted among the starbursts. 

In summary therefore, considering only those galaxies in 
which sufficient data is available, we have 12 galaxies with an 
unknown far-infrared excess relative to our models, 12 with 

AGN fractions /agn > 0.3, 16 starburst-dominated(/AGN < 0.3). 
Note that we previously found that AGN bolometric fractions 
above 0.3 cannot be reliably measured in this broad-band fit- 
ting (Negrello et al. 2009). If we exclude the poorly-fit SEDs, 
we find that (43 + 9)% of our submm-selected galaxies have 
/agn > 0.3; alternatively, if we conservatively include the galax- 
ies with far-infrared excesses as starburst-dominated, we find 
(30 + 7)% have AGN bolometric contributions above this thresh- 
old. We therefore treat the Icr range to be (23 - 52)%. Pope 
et al. (2008) used mid-infrared spectroscopy of submm-selected 
galaxies and found that a fraction (15 ± 10)% (i.e. 2/13 galaxies 
in their sample) had an active nucleus which contributes signif- 
icantly to ffie mid-infrared. This is slightly lower than found in 
our study, though only marginally inconsistent given the small 
number statistics in both studies. 

Eight of our SHADES galaxies have been identified as hav- 
ing anomalous 610 - 1400 MHz spectral indices (Ibar et al. 
in preparation): LOCK850.001, 015, 018, 022, 024, 033, 040, 
087. Of these, LOCK850.001, 024 and 033 have flat radio spec- 
tra (-dlnS'v/dln V < 0.5), while the remainder have unusually 
steep spectra, consistent with synchrotron ageing of electrons 
in the radio lobes from active nuclei. We have identified three 
of these eight securely as having AGN bolometric contributions 
/agn > 0.3, and three with /agn < 0.3. This suggests that ev- 
idence of AGN in the radio is not necessarily an indicator of 
bolometrically-dominant active nuclei. 

The high star formation rates (several hundred Mo/year, e.g. 
Hughes et al. 1998) and low number density of submm-selected 
galaxies (e.g. Scott et al. 2006, Coppin et al. 2006) suggest that 
submm galaxies are short-lived. Indeed the steep number counts 
of submm-selected galaxies necessarily imples that we sam- 
ple submm galaxies at around their peak phase of far-infrared 
luminosity. Even a slowly-varying monochromatic luminosity 
will induce this effect: in the models of Takagi et al. (2003, 
2004), the submm luminosity varies by around a factor of 2 - 3 
over a timescale of ~ 200 Myr. Denoting this for simplicity as 
Syit) oc e-0-5(^/<^)\ where T measures the secular evolution and 
cr - (200/2.35) Myr is the timescale of the variation, we find 
that a power-law differential source counts of dN/dS oc 5 " 
gives an observed age distribution in a flux-limited sample of 

Pr(T) cc e-("+3)(7)' = e-^-'^^y (4) 

where (cr')^ = cr^ /(2a+l). This is a far shorter timescale than the 
intrinsic variation, given the observed counts slope of a ^ 3.5-4 
(Coppin et al. 2006). 

There have been suggestions from numerical simulations 
(e.g. di Matteo et al. 2005) that feedback from active nuclei is ca- 
pable of truncating star formation activity. However quasars also 
have star-forming hosts (eg Serjeant & Hatziminaoglou 2009) 
which may pose challenges for models in which star formation 
is truncated too abruptly by active nuclei (e.g. Narayanan et al. 
2009). Therefore we suggest that the energy input from active 
nuclei through quasar-mode or radio-mode feedback does not 
immediately truncate star formation, but rather suppresses it on 
the same timescales as the quasar lifetime itself. We may also 
find that star-forming far-infrared-luminous populations selected 
at shorter wavelengths than SHADES have warmer colour tem- 
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peratures (e.g. Blain et al. 2003) and higher AGN bolometric 
contributions. In this interpretation, these would be later phases 
in the co-evolution of active nucleus and starburst. This would 
also be consistent with the observation that the /T-band Hubble 
diagram of hyperluminous starbursts is tight when the hyperlu- 
minous galaxies are selected at 60 /urn (Serjeant et al. 2003b), but 
has a much higher dispersion when the hyperluminous galaxies 
are selected at submm-wavelengths (e.g. Smail et al. 2004). 

Finally, our models predict submm-selected galaxies have 
fluxes of ^ 1 - 10 mJy at an observed frame wavelength of 70 fim. 
This was beyond the capabilities of both AKARI and Spitzer 
for direct detections, but is consistent with the Spitzer stack- 
ing analysis detection by Dye et al. (2007) of (3.63 + 0.77) mJy. 
Direct detections at this depth may be accessible to Herschel at 
shorter wavelengths, and would probe the hot dust components 
in submm-selected galaxies, currently not strongly constrained. 
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Fig. 6. SED fits to submm-selected SHADES galaxies in the Lockman Hole East, using models from Takagi et al. (2003, 2004) and 
Efstathiou & Rowan-Robinson (1995). 
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Fig. 7. SED fits to submm-selected SHADES galaxies in the Lockman Hole East, using models from Takagi et al, (2003, 2004) and 
Efstathiou & Rowan-Robinson (1995), 
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Fig. 8. SED fits to submm-selected SHADES galaxies in the Lockman Hole East, using models from Takagi et al, (2003, 2004) and 
Efstathiou & Rowan-Robinson (1995), 
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Fig. 9. SED fits to submm-selected SHADES galaxies in the Lockman Hole East, using models from Takagi et al. (2003, 2004) and 
Efstathiou & Rowan-Robinson (1995). 
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Fig. 10. SED fits to submm-selected SHADES galaxies in the Lockman Hole East, using models from Takagi et al, (2003, 2004) 
and Efstathiou & Rowan-Robinson (1995), 



